Reconstructing Baryon Oscillations : A Lagrangian Theory Perspective 
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Recently Eisenstein and collaborators introduced a method to 'reconstruct' the linear power spectrum from 
a non-linearly evolved galaxy distribution in order to improve precision in measurements of baryon acoustic 
oscillations. We reformulate this method within the Lagrangian picture of structure formation, to better under- 
stand what such a method does, and what the resulting power spectra are. We show that reconstruction does 
not reproduce the linear density field, at second order. We however show that it does reduce the damping of 
the oscillations due to non-linear structure formation, explaining the improvements seen in simulations. Our 
results suggest that the reconstructed power spectrum is potentially better modeled as the sum of three different 
power spectra, each dominating over different wavelength ranges and with different non-linear damping terms. 
Finally, we also show that reconstruction reduces the mode-coupling term in the power spectrum, explaining 
why mis-calibrations of the acoustic scale are reduced when one considers the reconstructed power spectrum. 

PACS numbers: 



I. INTRODUCTION 

The baryon acoustic oscillation (BAO) method 01] is an in- 
tegral part of current and next-generation dark energy experi- 
ments. Oscillations in the baryon-photon fluid, frozen into the 
matter distribution at decoupling, provide a standard ruler to 
constrain the expansion of the Universe. These sound waves 
imprint an almost harmonic series of peaks in the power spec- 
trum P(k), corresponding to a feature in the correlation func- 
tion £(r) at ~100 Mpc, with width ~ 10% due to Silk damp- 
ing (see HI for a detailed description of the physics, and 
Jj] for a comparison of Fourier and configuration space pic- 
tures). While the early Universe physics is linear and well 
understood, the low redshift observations are complicated by 
the non-linear evolution of matter (not to mention galaxy bias 
and redshift space distortions [5], but we will defer these to 
future work) whicherases the oscillations on small scales and 
shifts the peaks 



P ohs (k) = e~ k2 * 2 / 2 PUk) + Pmc(k) + 



(1) 



by coupling individual fc-modes which are at early times inde- 
pendent. The exponential damping of the linear power spec- 
trum (or equivalently the smoothing of the correlation func- 
tion) reduces the contrast of the feature and thereby the preci- 
sion with which the size of ruler may be measured. Neglect 
or incorrect modeling of the "mode-coupling" term P mc may 
bias the resulting distance measurements. 

In 10] it was pointed out that much of the modification 
to the power spectrum comes from large-scale modes, bulk 
flows and super-cluster formation, in principle enabling their 
effects to be corrected. Eisenstein et al 19[] introduced a 
method for removing the non-linear degradation of the acous- 
tic signature, sharpening the feature in configuration space 
or restoring/correcting the higher k oscillations in Fourier 
space. Given the ambitious nature of future experiments, there 
has been considerable interest H [l(| [Htl in "reconstruction" 
schemes which remove the effects of non-linearities, reducing 
the damping and mode coupling terms above. 



Since the method proposed in \$\ is an inherently non-linear 
mapping of the observed density field, it is difficult to intu- 
itively understand. It is however easily formulated within the 
Lagrangian picture of structure formation, where the funda- 
mental quantity is the displacement of particles from their ini- 
tial positions (contrasted with the Eulerian picture where one 
tracks the evolution of the density field at a fixed location). 
Motivated by recent developments in Lagrangian perturbation 
theory (LPT) J3, fl2h . we discuss reconstruction within the 
context of LPT, both to elucidate how it works and to expose 
possible shortcomings. Although we use the method of 131 for 
specificity, the lessons learned have broader validity. 

We proceed as follows : ^introduces the essential aspects 
of both LPT as well as reconstruction. We then compute the 
reconstructed density field to second order, and demonstrate 
that there are corrections to the linear density at this order. ^UJ\ 
then explains why the BAO feature is enhanced in the recon- 
structed power spectrum. We conclude in £|IV1 highlighting 
potential avenues for improvements. 



II. RECONSTRUCTION AND THE DENSITY FIELD 

The Lagrangian description of structure formation |[l3l [Til 
[l5tl relates the current, or Eulerian, position of a mass ele- 
ment, x, to its initial, or Lagrangian, position, q, through a 
displacement vector field ^(q), 

x = q+*(q). (2) 

The displacements can be related to overdensities by flrjll 

(5(x) = J d 3 q (5 (D) (x q *) 1 . (3) 

where is the 3D Dirac 5 function, or in Fourier space by 

5(k) = J d 3 q e" ik q (V lk * (q) - l) . (4) 
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The displacements evolve according to 

[q+*(q)] 



d 2 * „ 



(5) 



where is the gravitational potential. Analogous to Eulerian 
perturbation theory, LPT expands the displacement in powers 
of the linear density field, Si, 



iff = VJ/W -+- *J>( 2 ) 



(6) 



where 111 

«r(»)(k) = 




x LW(k 1 ,...,k n ,k)5 I (ki)---Mkn 



(7) 



and the LA") have closed form expressions, generated by re- 
currence relations. Specifically, 



(8) 



is the well known Zel'dovich displacement [e.g. [181], which is 
1 st order LPT. Expanding the exponential in Eq. we obtain 
a perturbative series for the overdensity, <5 = S 1 - 1 ^ + S^ + ••• 
where, e.g., 



5^ = I d 3 q 
or in terms of the L^ n ''s 



ik* (2) 



(k.*W) 2 



(9) 



^|^ (D)(ki+k2 _ k) 

x 5i(ki)Si()L2) [k-L( 2 )(k l7 k 2 ,k) 
+ k-LW(ki)k-LW(k 2 ) 



(10) 



is second order in the linear density field Si. 

The prescription of |§] can be cast into this framework as 
follows: 

• Smooth the density field to filter out high k non- 
linearities. In Fourier space, this is equivalent to mul- 
tiplying by a function S(k) which monotonically de- 
creases from unity at low k to zero at high k, 



6Ql) -> S(k)S(k) . 



(11) 



• Compute the negative Zel'dovich displacement from 
the smoothed density field 



(12) 



• Shift the original particles by s and compute the "dis- 
placed" density field, 

S d (k) = y> ge -* k <l (V*M*(q)+s(q)] _{j . (13) 

Note that if the original density field were linear, and 
5 = 1, this would undo their displacements exactly, 
moving them back to their original positions and giving 
S d = 0. 

• Shift an spatially uniform grid of particles by s to form 
the "shifted" density field, 



S.Qt) = I d 3 qe- ik * ( e - *'"^ - 1 



(14) 



Again, assuming linear theory would imply <5 S (k) 
-<J(k). 

• The reconstructed density field is defined as S Tecon 
Sd - S s 

^rccon(k) = J d 3 9 e-* k q e- <k S ( q ) (e- 4k * (q) - 1 

with power spectrum P rccon (fc) oc (|<5 2 ec0 n|)- 



(15) 



Note that S oc s — > is equivalent to no reconstruction, which 
is helpful in intepreting some of the expressions below. 

When applied to simulations this process yields an en- 
hanced BAO feature [8l |9j [lOt \lM with a reduced 'shift' in 
the peak. Our focus here is to understand what this procedure 
is doing within an analytic framework. 

In the spirit of LPT we can expand the reconstructed density 
field in a perturbative series 



5 - 

Brecon — "rccon 



<J(2) 
rccon 



(16) 



As anticipated above, the reconstructed field equals the linear 
density field to lowest order. Working to the next order, we 
find 



<J(2) 
w rccon 



5 (2) 



2 J (2tt) 3 
x <5 / (k 1 ),5 / (k 2 )k-L( 1 )(k 1 )k-L( 1 )(k 2 ) 
x [5(kx)+5(k a )] . (17) 



We observe that the second-order term in the reconstructed 
density field does not vanish. While S^ contains L^ 2 ), the 
correction only involves and so cannot fully cancel the 
non-linearity. This is a general feature - the corrections to 
on iy involve terms \J- l<n ^> - and follows from the fact 
that we only worked to first order when shifting objects. We 
note in passing that one might be able to construct higher or- 
der reconstruction schemes such that (5 ; ( '™ >1 ^ contributions to 
the reconstructed density vanish, but that is beyond our scope 
here. 

To recap: the reconstruction algorithm above generates a 
density field with second order corrections, not the linear den- 
sity field. The next section explains why simulations saw an 
improvement when using reconstruction, by considering the 
reconstructed power spectrum. 
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III. THE POWER SPECTRUM 



A. A Toy Model 

To best highlight the effects of reconstruction on the power 
spectrum, we start with a toy model that captures both the 
physics and the algebraic structure of the full gravitational 
perturbation problem. This toy model is particularly useful 
for identifying the effect of reconstruction on the nonlinear 
damping of the linear power spectrum in Eq. Q] § MI B I de- 
scribes the correspondence between the toy model and the full 
gravitational instability problem, extending the analysis of the 
effect of reconstruction to the mode coupling terms as well. 

Consider a model, inspired by the peak-background split, 
where \^ can be split into low (L) and high {H) frequency 
pieces, 



10 



o 



w 



(18) 




0.01 



10 



(h/Mpc) 



with \P l the Zel'dovich displacement based on a linear den- 
sity field Si, 



*i(k) 



<£*«(k) 



(19) 



For simplicity we assume that is also Gaussian and is 
uncorrected with l . The intuitive picture behind this model 
is that l encodes the linear density field, while encodes 
the non-linearities; importantly, the baryon oscillations only 
exist in \I>l and not in ^ h ■ 

Using Eq. (HJi the power spectrum is 



P(k) = J , 



qe 



- zkq 



-it,A*,(q) 



- 1 



(20) 



where q = qi — q2, and A\& = ^(qi) — , 4'(q2). For Gaus- 
sian 4' 



e " k A *^ ) = exp 



■-kikj (A* i (q)A%(q)} 



(21) 



with 



hk 3 (A^(q)A^(q)) = 2fc 2 <* 2 (0)) - 2fc i *,-&(q) , 

(22) 

where £ij(q) = iiox)^ is the displacement correla- 

tion function and we have used translational invariance for the 
correlation function at zero lag. To lowest order the zero-lag 
correlation function is £y(0) = <5jj£ 2 /2, with £ 2 the mean- 
squared Zel'dovich displacement of particles, 



^JdpP L (p) (23) 



with a similar expression for . Note that the relation of 
the damping to the Zel'dovich displacement follows naturally 
from the LPT formalism and shows the similarity of the treat- 
ments in Refs. fiSd and fl. 

Given our assumption of uncorrelated low and high fre- 
quency pieces, we have S 2 = S| + Y? H . However Fig. Q] 



FIG. 1 : The damping scale at z — as a function of the maximum 
wavenumber for the linear (solid) and non-linear (dashed) power 
spectra. Note that the dominant contribution to the damping scale 
come from linear motions. 



demonstrates that the dominant contribution comes from rel- 
atively large (k < 0.3/iMpc" 1 ) scales. If we estimate 
by substituting the non-linear power spectrum in the equation 
above, we find that the dominant contribution comes from lin- 
ear motions, even at z = 0. For simplicity, we will therefore 
assume £ 2 ~ Y? L in what follows. 

The non-linear power spectrum is then given by 



P(k) 



-k 2 T. 2 L /2 



d 3 q e" 



(24) 



Following 0] we leave the zero-lag piece exponentiated, but 
expand the exponential inside the integral. The first term of 
the expansion gives P^. This procedure can be viewed as a 
resummation of terms in the standard perturbative expansion 
which leads to a power spectrum of the form in Eq. ((T), 



,Pob s (fc) = e- k "^ /2 P L (k) + P mc (k) + 



(25) 



with Pi in = P L (k). Note that P mc {k) contains terms 0{^> 2 H ) 
representing the high frequency part of the power spectrum 
and terms 0{^f\) corresponding to second order (in Pl) cor- 
rections. We will consider these terms in the next section. 

The above can be extended to compute the reconstructed 
power spectrum for this model. Since * h has no low fre- 
quency piece by construction, we assume that the inferred 
shift, s(k), is simply given by 



s(k) = -S(fc)* i (k) + 0(* 2 ) 



(26) 



The fields 5d and S a of Sec. [n]are then generated to first or- 
der by (1 — S)^l + and — S^l respectively. Since 
the reconstructed density field is the difference of the two 
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fields, there are three terms (two auto-spectra, P ss ,Pdd and 
one cross-spectrum P s( i) that make up the reconstructed power 
spectrum: P r ocon = Pss+Pdd—^P s d- The auto-power spectra 
are exactly analogous to the non-linear power spectra, except 
for the damping terms, 

P ss (k) = e- k2 ^/ 2 S 2 (k)P L (k) + ■■■ (27) 

and 

Pdd(k) = e- fe2s -/ 2 [1 - S(k)} 2 P L {k) + ■■■ (28) 

where we've dropped higher order terms. The Gaussian 
damping is modified to 

S - = 3^2 / d P S 2 (p)P L (p) (29) 

with an analogous expression for Edd with S 2 — > (1 — S) 2 . 
The cross power spectrum is 

P sd (k) = ~e~ k2 ^l 2 S(k)[l - S(k)]P L (k) + ■■■ (30) 
where 

sL = ^(sL + sL) , (3D 

and the negative sign comes from the fact that the random 
field was shifted by the negative Zel'dovich term. Putting the 
pieces together, we find that the damping term becomes 

D(k) = e - fc2s2 /2 _> 5 2 (fc) e - fe2s "/ 2 

+ [1 - S{k)] 2 e- k ' 2 ^' 2 

+ 25(fc)[l - S(k)}e- k2 ^ /2 .(32) 

Before proceeding, it is useful to choose an explicit form 
for the smoothing; the standard choice is a Gaussian, 

S(k)=exp( . (33) 

Fig. |2] plots the various damping scales as a function of the 
smoothing scale. As expected, for non-zero smoothing, both 
S ss and Y.dd (and therefore T, s d as well) are less than the non- 
linear damping scale. This is the crux of the reconstruction 
method - that the Pl contribution to the reconstructed power 
spectrum is less damped than in the nonlinear power spec- 
trum. This holds even when taking into account that there 
are additional terms depending upon S(k) in D(k) as we now 
show. 

Before considering Eq.[32]for arbitrary choices of smooth- 
ing scales, we consider the special case where S ss = T,dd = 
S s d; for the Gaussian smoothing above, this corresponds to 
a smoothing scale R ~ 30 /i _1 Mpc. The damping of Pl in 
the reconstructed power spectrum simplifies considerably; the 
reconstructed power spectrum has the form, 

Prccon(fc) = e- fe2E -/ 2 P L (fc) + • • • (34) 




10 20 30 40 50 
R (Mpc/h) 



FIG. 2: The ratio of E ss , E s( ; and E<m to E, as a function of the 
Gaussian smoothing scale, R. Note that for no smoothing, E sa = E 
and E^ci = 0, while for infinite smoothing, E^d = E with E ss = 0. 

Note that this is identical to the form of the nonlinear power 
spectrum (Eq. l25l l except that E ss < S, reducing the damp- 
ing. 

Fig.[3]shows the damping, D(k), for smoothing scale R = 
5 /i _1 Mpc, as an example of its general form for an arbitrary 
choice of smoothing scale. Given R, the factors involving 
S(k) determine the range of wavenumbers for which each of 
the three power spectra dominate. For large R (see below), 
Pdd dominates over the wavenumbers important for baryon 
oscillations (0.07 < fc(/i/Mpe) < 0.35) but as we argue be- 
low, this limit is not optimal. As we decrease R, we might 
have expected that P ss would have dominated; however, de- 
creasing R quickly increases E ss to close to the nonlinear 
damping scale (Fig. [2j>, limiting the importance of P ss . In- 
deed, in Fig. [5J we see that P ss has the linear power spec- 
trum more strongly damped than the nonlinear power spec- 
trum. The dominant term at small R is therefore P s d\ Fig. [2] 
shows that Y, s d is ~ 0.6 S and is only weakly dependent on 
R. This suggests that such a reconstruction method can reduce 
the damping of the linear power spectrum by a factor ~ 2. 

The above discussion argues that the smoothing scale de- 
termines the wavenumber where Pdd becomes dominant. The 
obvious question is whether the above analysis suggests a 
value for the smoothing scale. We argue that the natural 
choice is R ~ S, the nonlinear (damping) scale. To see why, 
we start by observing that the terms we ignored in Pdd are 
0(^ 2 H ), whereas for P ss and Pdd they involve higher powers 
of the displacement. This is just the statement that the small- 
scale displacements have their largest effect on Pdd which is 
not surprising, given that Pdd is based on the original density 
field. We would ideally want to reduce these terms, which 
argues for making R as small as possible. However, from 
Eq.|25] we see that the linear field is damped on scales smaller 
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FIG. 3: The damping of the linear power spectrum for the nonlinear 
power spectrum (dashed line), and the reconstructed power spectrum 
(Eq.f32] solid line, assuming a smoothing scale R = 5/i _1 Mpc). 
The dotted lines decompose the reconstructed damping into the lead- 
ing contributions from its P sa , P s d and Pdd components These 
curves have been calculated assuming z = 0. 



than S. Smoothing on scales much smaller would then violate 
our assumption that s(k) is derived from the linear density 
field, which leads to choosing R ~ S as might have intu- 
itively been expected. 

The above discussion explains how reconstruction reduces 
the damping of the acoustic oscillations (or equivalently, how 
it sharpens the peak in the correlation function). We now turn 
to its effect on the mode-coupling terms, by considering the 
reconstructed power spectrum within LPT. 



B. Lagrangian Perturbation Theory 

Many of the features of reconstruction in the last section 
carry across to the gravitational instability problem within 
LPT. We will closely follow the LPT formalism developed 
in (|7l [HI] in which the broadening of the peak and the mode 
coupling terms appear naturally. 

For the un-reconstructed power spectrum the derivation 
leading to Eq. (l20b still holds. However now we must use 
the cumulant expansion theorem 



<e-«) 



exp 



E 

N=l 



H) 

TV! 



N 



(35) 



(where the (X N ) C are the connected moments) to compute 
the expectation value of the exponential. In the toy model 
only the N = 2 term survived, for the full problem higher 
orders contribute as well. Expanding (k ■ A\E') Ar using the 
binomial theorem we have two types of terms: those where 



the * are all evaluated at the same point (which we can take 
to be the origin) and those with iqis and N — j q 2 s. As 
in the toy model, and following 0], we leave the first set of 
terms exponentiated while expanding the second set of terms 
in powers of \&. If we keep only the lowest order terms in the 
exponential we regain the form of Eq. (Q~|) with £ given by the 
rms Zel'dovich displacement 



P(k) 



3 -fc 2 S 2 /2 



P L {k) 



1 + j rf 3 fc 1 P L (fc 1 )G(k,k 1 ) 
d 3 k 1( i 3 k 2 P L (k 1 )P i (k 2 ) J F 2 (k 1 ,k 2 ,k) 



(36) 



where F and G can be expressed in terms of L^ 1 ) and L^ 2 ' 
and explicit expressions may be found in 10]. There are os- 
cillations in Pl and the mode-coupling term (third line), but 
the integral in the second line has a wide kernel so the os- 
cillations are suppressed. As it happens, F is peaked around 
ki w k 2 « k/2, which helps to explain why this the term in 
the third line leads to a peak shift. If Pi contains an oscilla- 
tory piece, e.g. sin(fcr), then the third term contains a piece 
schematically of the form sin 2 (fcr/2) ~ 1 + cos(fcr), which 
oscillates out of phase with Pl. It is the sum of the two out of 
phase components that leads to a shift in the peak of £ (r) or 
the phasing of the harmonics in P b s (k)- 

It is now straightforward, though tedious, to repeat these 
steps for the reconstructed field. The formalism of Ref. J3l 
must be generalized to allow two displacements (s and \E'). 
Again there are three contributions, P ss , Pdd and P s d, and 
three smoothings, S ss , Y,dd and £ s d of the same form as be- 
fore. The term proportional to Pl becomes 



= {e- fe2s -/ 2 S 2 (fc) 

+2e~ k2 ^/ 2 S(k) [1 - 5(A)] 

+ e -fe 2 sL/2 [l-5(fc)] 2 }p L (fc). 



(37) 



directly comparable to the result of the toy model. 

The leading contribution to the mode coupling term is the 
same as in standard perturbation theory, and is strictly posi- 
tive, coming from (S^S^y Recalling the relation between 

Brecon and we need to replace F in the mode-coupling 
term with F s.t. 

2P(k 1 ,k 2 ) = k-L( 2 )(k 1 ,k 2 )+k.L( 1 )(k 1 )k-L( 1 )(k 2 ) 
x [l-5(k!)-5(k 2 )] (38) 

where ki + k 2 = k. The piece of the mode-coupling integral 
which shifts the peak comes from ki w k 2 w k/2 where 
k • L 2 — and 



P(k/2,k/2) = 2{l-25(k/2)} . 



(39) 



Since the term in {• • •} is bounded between —1 and 1, F 2 < 
F 2 for all k, suppressing the reconstructed mode-coupling 
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term relative to the corresponding term in the non-linear 
power spectrum. This explains why the mis-calibrations in 
the acoustic scale were reduced after reconstruction in iH . 

IV. COMMENTS 

It is now generally understood (0>H01 an d this work) that 
the dominant effect of the non-linear evolution of matter per- 
turbations on the baryon oscillations is to damp the higher har- 
monics, Pobs(fc) = exp(— fc 2 E 2 /2)Pi; n (fc) + • • -, or equiva- 
lently, smooth the feature in the correlation function. Eisen- 
stein et al U proposed a "reconstruction" method, demon- 
strated on simulations, that undoes this non-linear smoothing 
and appears to restore the linear power spectrum. Motivated 
by recent progress in Lagrangian perturbation theory lffl [l2tl . 
we revisit this algorithm in order to better understand why it 
works as well as its shortcomings. Our principal conclusions 
are 

(i) The field generated by the reconstruction process is not 
the linear density field at second order. Note that this is 
a general statement, independent of assumptions about 
the smoothing of the initial density field. 

(ii) Reconstruction does reduce the damping of the oscilla- 
tions, by about a factor of 2 when the input density field 
is smoothed on the non-linear scale. 

(iii) Reconstruction also reduces the mode coupling terms 
which introduce an out of phase component of the os- 
cillations or shift the peak. 

(iv) The reconstructed power spectrum is the sum of three 
power spectra (the auto-power spectra of the displaced 



and shifted fields, and their cross-spectrum), each of 
which have different damping terms (Eq. [32l . An ap- 
propriate model for the reconstructed power spectrum 
should take this into account, instead of modeling it as 
a single damping scale. 

(v) When the smoothing scale is close to the non-linear 
scale, the correlation between the shifted and displaced 
fields plays a crucial role. 

Our results suggest a number of natural extensions. The 
effects of bias and redshift space distortions have been in- 
corporated into the Lagrangian formalism fK. [l2Tl . and could 
therefore be folded in to the LPT formulation of reconstruc- 
tion. We have observed that the reconstructed density field 
is not the linear density field; an interesting possibility is to 
explore whether higher order reconstruction schemes actually 
yield dividends. Even within the context of the existing recon- 
struction schemes, it is possible that a different weighting of 
the three power spectra may yield improved accuracy in mea- 
suring the distance scale. We leave these avenues open for 
future investigation. 
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